Extended wave propagators as pulsed-beam 
communication channels 

Plenary lecture at Days on Diffraction Conference, St. Petersburg, Russia 
May 30-June 2, 2006. http://math.nw.ru/DD/ 

Gerald Kaiser 

Physics Department, UT- Austin 
and Center for Signals & Waves 
kaiser@wavelets.com 
http : //wavelets . com 

February 7, 2008 



Abstract 

Let P{xr — Xe) be the causal propagator for the wave equation, 
representing the signal received at the spacetime point Xr due to an 
impulse emitted at the spacetime point Xe- Such processes are highly 
idealized since no signal can be emitted or received at a precise point in 
space and at a precise time. We propose a simple and compact model 
for extended emitters and receivers by continuing P to an analytic func- 
tion P{zr — Zf,), where Zf. = Xe + iye represents a circular pulsed-beam 
emitting antenna centered at Xe and radiating in the spatial direction 
of j/e while Zr = Xr — iVr represents a circular pulsed-beam receiv- 
ing antenna centered at Xr and receiving from the spatial direction of 
yr. The space components j/ei l/r of !/e, yr give the spatial orientations 
and radii of the antennas, while their time components Se, Sr represent 
the time a signal takes to propagate along the antennas between the 
center and the boundary. The analytic propagator P represents the 
transmission amplitude, forming a communication channel. Causality 
requires that the extension/orientation 4- vectors and yr belong to 
the future cone V^, so that Zf, and Zr belong to the future tube and the 
past tube in complex spacetime, respectively. The imaginary "retarded 
time" Te = — lUel/c represents the duration of the emitted pulse, 
and Tr = Sr — ll/rl/c represents the integration time for the received 
pulse. The bandwidths of the antennas are l/Tg and l/T^, respec- 
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tivcly. The invariancc of P under imaginary spaeetime translations 
(ze — > Zg + ir], Zr —>■ Zr + irj) has nontrivial consequences. 

1 Introduction 

The purpose of this paper is to consider the following proposition: When 
extended analytically to complex spaeetime in a sense to he made precise, 
retarded wave propagators describe the causal interaction between a source 
of finite spatial extent and a sink of finite spatial extent. The supports of 
the extended source and sink are give by the imaginary spaeetime coordi- 
nates associated with the analytic continuation, up to an equivalence due to 
invariance under complex spaeetime translations. 

We shall explain and prove this proposition in the case of the scalar wave 
equation. It also extends to Maxwell's equations [6, 7, 8] and possibly some 
other systems [2]. Begin with the retarded wave propagator 

6(t — r) 

pi^x) = — -, X = Xr - Xe = {X,t) eW^, r =\x\, (1) 

Airr 

where we have set the constant propagation speed to c = 1 by an appropriate 
choice of units. P is a fundamental solution of the wave equation, 

DP{x) = {df - A)P{x,t) = d{x)6{t) = 6{x). 

If we fix the emission event Xe = {xe,te) and vary the reception event 
Xr = {xr,tr), then P{xr — x^) may be regarded as a wave emanating from 
the emission point Xe at time te, the wave front at any time tr being the 
sphere \xr — Xe\ = tr — te (with no wave when tr < te)- This gives the usual 
picture of an exploding wave propagating causally away from x^- 

However, we can also fix the reception event Xr and examine P{xr — x^) as 
a function of the emission event Xe- Then the same causal propagator is 
seen to describe imploding wave fronts \xe — i c^l — — {te — ^7^) ^t emission 
times te < tr, absorbed at Xr at te = tr with no waves at te > tr- Although 
this situation is completely symmetrical to the one with fixed emitter, it 
causes a bit of discomfort because we never actually see waves imploding 
toward a point. But a little reflection shows that we do not "see" acoustic 
or electromagnetic waves exploding either (like waves on a pond). Rather, 
the behavior of P{xr — Xe) as a function of Xr is measured by setting up a 
system of idealized "test receivers" (a spaeetime version of point-receivers) 
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at various x,. and piecing together their outputs. When viewed hke this, 
the imploding wave picture no longer seems in conflict with causality: fix a 
single point receiver at Xj. and set up a system of "test emitters" at various 
Xe- The resulting transmission is again P(xr — Xe), this time as a function 
of Xe- Note that applying of the wave operators Dg, Dr in the emission and 
reception variables Xe,Xr, respectively, gives 

\I\eP{Xr — Xe) = \I\rP{Xr — Xg) = 5{x), P{Xr — Xe) = for > tr, 

confirming that P{xr — Xe) is indeed an advanced wave propagator in Xe- 

The reason why fixing the receiver and varying the emitter seems unusual 
is that the field concept itself is biased toward reception, as proved by its 
very terminology: the filed variable is the "observation point," never the 
"emission point." The bias toward reception also manifests itself in the fact 
that while no one has a problem with emission from a extended source, 
the dual process of observing with an extended receiver is usually treated 
indirectly by resorting to abstract reciprocity theorems involving two anti- 
causal twists: first the emitter and receiver are exchanged, then a space-time 
reversal is performed (represented by a complex conjugation in the frequency 
domain). By contrast, the above version of reciprocity (focus on Xg instead 
of Xr) is direct and never leaves the causal world. It automatically involves 
a space-time reversal since x = x^ — Xg. 

We shall treat extended emitters and receivers in a symmetric manner by 
showing that an analytic extension of P{xr — Xe) to complex Xg and Xr 
represents emission and reception by extended objects that, in fact, resemble 
antenna dishes whose sizes and orientations are given by the values of the 
imaginary spacetime variables attached to Xg and Xr, respectively. 

2 The extended wave propagator 

There is, of course, no way to analytically continue the numerator 6{t — r) 
of (1). But it follows from the Sochocki-Plemelj formula that 



hence P{x, t) is the jump in boundary values of the Cauchy kernel 





6{t - is) 
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(3) 
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going from the lower complex time half-plane to the upper half-plane. Thus 
we begin by extending P{x) to complex time, leaving space real for now: 

6(T-r) 1 
P(x,t) = —^ - = ^— -, T = t-is. (4) 

Next, we complexify space by replacing the Euclidean distance r = \x\ by 
the complex distance [4] 

f{z) = \^ = y/z ■ z, z = X — iy ^ 



= \/r2 - a2 - 2mx3, a=\y\, X3 = xy, (5) 

where y = y/a and it is assumed that a > 0. f is an analytic but double- 
valued function in C'^ which will be restricted to a double- valued function in 
by fixing y 7^ 0. Substituting r ^ f in (4), we obtain the full extension 
of P to complex spacetime: 

P{z) = P{z, t) = ^^^^ = -, z = (z, t) G C\ (6) 

This will be the main object of our attention. We will show that while the 
boundary values of P in the sense of (2) describe idealized emission and 
reception events, the interior values of P{zr — Ze) describe more realistic 
processes representing emission from and reception in extended spacetime 
regions surrounding the events Xe = Re Ze and Xr = Re z^. The imaginary 
spacetime variables ye = Im Ze and yr ^ — Im Zr will be seen to determine 
the extent of the spacetime regions where the emission and reception occur. 
They effectively give the scales (durations and source dimensions) in space 
and time for the emission and reception processes, subject to equivalence 
under complex spacetime translations z'^ = z^ + C: z'^. = Zr + C- 



3 Interpretation of P as a pulsed beam 

To interpret the extended propagator P, we begin by studying the complex 
distance f. Writing 

r = p{x) - iq{x), (7) 

where the dependence on the fixed imaginary variables y = — Im z is sup- 
pressed, so (5) gives 

p —q =r —a and pq = ax^. (8) 
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Hence the cylindrical coordinate p orthogonal to y is given by 

2 2 2 2,2 2 22/2 {o? + p'^){o?' — q'^) , - 

p =r -xn = a +p -q -pq/a = ^ . (9) 

This shows that q^ < , with equality only on the 2:3 axis, and 

2 2 / 2 2\ ^ 2 
p = r — [a — q ) <r , 

with equality again only on the X3 axis. Thus f satisfies the bounds 

\p\ = I Re f{z)\ < I Re z\ = r and \q\ = \ Im r(z)| < | Im z| = a, 

with equality if and only if x is parallel to y. Furthermore, (9) and (8) give 

+d = l P' -d = i (10) 

showing that the level surfaces of p'^ are oblate spheroids Sp, and those q"^ 
are one-sheeted hyperboloids TCq. These are easily proved to be orthogonal 
families sharing the circle 

C = {x r = a, xs = 0} (11) 

as a common focal set. Since C is precisely the set in M'^ where f = 0, 
it is also the set of branch points of r in (for given y ^ 0). When f 
is continued analytically within from any point x ^ C around a simple 
closed loop threading C, it changes sign due to the double- valuedness of the 
complex square root. To make r a single- valued function in M^, we choose 
a branch cut B consisting of a membrane spanning C, 

dB = C, 

and declare that r changes sign upon crossing B, so that continuation along 
a closed loop now returns f to its original value. The simplest choice of 
branch cut is the disk 

£o = {xe M^l p = 0} = {:e G M^l r < a, = 0}, 

which is the degenerate limit of £p as p ^ 0. We choose the branch of f 
with p > 0, so that in the far zone (5) gives 

r ^> a f^r — iacosO, cos6 = x-y. 
5 



Hence (6) gives the far-zone expression 



r>a ^ ^ — (12) 

which is seen to be a pulsed beam of 0-dependent duration 

T{e) = s-acos9 (13) 
propagating in the direction oi y [6 = 0) whose peak value occurs at t = r: 

^P"^'' 87r2r(s-acose) " 8Tr^rT{e) ' ^^^^ 

To avoid singularities away from the source, we must require that the imag- 
inary spacetime coordinates satisfy the constraint 

s>a = \y\, i.e., y = {y, s) £ V+ = {{x,t) £ R'^\ \x\ < t}. (15) 

That is, y must belong to the future cone V+ of spacetime, or z = x — iy 
must belong to the past tube, also known as the forward tube in quantum 
field theory [10]. We postulate that s is the time required for a signal to 
propagate along the antenna between the center and the boundary, so that 
the condition s > a = a/c is simply a statement of causality. From (14) we 
find that the far-field radiation pattern is 



87r2(s — a cos ( 



which describes an ellipse of eccentricity a/s with one focal point at the 
origin. Thus J- peaks sharply around ^ = as s^a"*". Note that P has no 
sidelobes. 



4 General pulsed-beam wavelets and their sources 

The pulsed beam P decays slowly in time and angle (0(t~^) and 0{{s — 
acos9)~^), respectively) because these decays come from the Cauchy kernel 
(3). Faster decays in both time and angle are obtained by convolving P with 
distributional signals. Choosing a driving signal go{t), we define general 
pulsed-beam wavelet by 

Wiz, t) = r dt' go{t')P{z, T-t') = ^-^^fj^ (16) 
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where g is an analytic signal of go, defined by 

1 P go{t')dt' . 
9{r) = — — , T = t-is. (17) 

W is interpreted as the signal broadcast by the extended source using go 
as input, goit) should have reasonable decay so that the integral converges, 
but it need not be analytic. Then g(r) is analytic everywhere except for the 
support of go: 

drgir) =0 Vr ^ supp^o C M. 

In particular, g is analytic in the upper and lower half-planes, given there 
by the positive and negative frequency parts of go, respectively: 

git - is) = r dw e(^s)e— )5o(u;), (18) 

where is the Fourier transform of go and is the Heaviside function. (This 
is a special case of the multidimensional analytic-signal transform [2, 3, 6].) 
If go{t) vanishes on any open set, the analytic functions ^(t) in the lower and 
upper half-planes are analytic continuations of one another. It is possible to 
achieve excellent localization in both time and angle, still without sidelobes, 
by choosing an appropriate distribution for go- For example, 

goit) = d^6{t) g{T) = df6{T) - 



The pulsed-beam wavelet (16) reveals the mathematical nature of the exten- 
sion P ^ P. Fixing any vector y G the jump discontinuity of W{x — iy) 
across real spacetime is 

lim {W{x + iey) - W{x - iey)} - ~ 



*0 + 47rr 

The right side is a hyperfunction, a distribution in representing the jump 
in boundary values of a holomorphic function in a complex spacetime domain 
[1,9], in this case the jump from the past tube to the future tube. Physcially, 
it represents the wave caused by broadcasting the original driving signal from 
an idealized point source: 

So{x, t)^a = goimx). 
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Thus we define the extended source of W{z) by applying the real wave 
operator to W{x — iy): 



where the subscript is a reminder that we are differentiating only with re- 
spect to X since the imaginary source coordinates are constant. (They char- 
acterize the source's radius a, orientation y, and duration T(0) = s — a along 
the beam axis.) Note that W{x — iy) is singular along the branch circle C 
(where f = 0) and discontinuous across the branch cut B (where f changes 
sign). Hence the wave operator Ox must be applied in a distributional sense 
[4, 6]. A formal differentiation of W gives 5 = 0, which shows that S is 
actually a distribution supported on the singularities of W, which form the 
world-tube in spacetime swept out by the branch cut B C M^. The behavior 
in time is analytic because r — f is bounded away from the real axis as long 
as y £ V^. So the important fact for us is that at any time t, the spatial 
support of the source S is the branch cut B. The distributional source S and 
its regularized versions (equivalent Huygens sources supported on surfaces 
surrounding B) have been computed in [6], both in the spacetime and the 
Fourier (wavenumber-frequency) domains. 

The physical significance of the analytic signal g{T — r) is that an extended 
source cannot radiate an arbitrary signal with infinite accuracy because 
waves from different parts of the source interfere. The factor e~^'^ in (18) 
shows that g{t — is) can faithfully reproduce time variations only up to the 
time scale s, and similarly g{T — f) can reproduce time variations only up 
to 0{s — q). This gives a physical interpretation of W and shows that the 
imaginary spacetime vector y = (y, s) is indeed a multidimensional version 
of the scale parameter in ordinary wavelet analysis [3] . 

5 Pulsed-beam wavelet communication channels 

Our main goal here is to show that the pulsed-beam wavelets W{x — iy) 
form natural communication channels. For this purpose, write the complex 
spacetime variables as 

z = x — iy = Zr — Ze , where Ze = Xe + iye and Zr = Xr — iyr , (20) 



S{x - iy) = \JxW{x - iy) 



(19) 



with 



) G V+ and yr = (y,., Sr) G 



+ ) 
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that is, 



Se > ae= IVel and Sr > ttr = \yr\- 



(22) 



Again, these are merely statements of causality if we postulate that Sg is 
the time it takes a signal to propagate along the emitting antenna from the 
center to the boundary and Sr is the time it takes a signal to propagate 
along the receiving antenna from the boundary to the center. The opposite 
signs in Zg and Zr are mathematically necessary to ensure that 



belongs to the past tube (as it must) because V+ is convex. The physical 
interpretation of the sign change is that while extended emission sources 
are parameterized by the future tube Xe + iye, extended reception sources 
are parameterized by the past tube Xr — iyr- (Roughly, emission is to the 
future while reception is from the past!) Though we shall not justify the 
definitions (20) further here, we now explore their consequences and show 
that they make complete physical sense. We shall interpret W{zr — Ze) as a 
transmission amplitude for a pulsed-beam wavelet emitted by an extended 
source centered at Xg and characterized by ye, then received by an extended 
source centered at Xr and characterized by y^- A given pair of source pa- 
rameters Zg and Zr thus represent the communication channel depicted in 
Figure 1. We have sketched the emitting branch cut Be, centered at Xe with 



radius Cg and orientation ye, and the receiving branch cut Br, centered at 
Xr with radius a,, and orientation y^. Note that the vector points out 
of Be while y,, points into Br- This is a graphical explanation of the sign 
difference in (20). To show that Figure 1 makes physical sense, note from 



Z = {xr - iyr) - {Xe + iye) = {Xr " Xe) - i{yr + ?/e) 




Figure 1: W{zr — Ze) communication channel. 
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(18) that g decays with increasing distance from the real axis. Hence we 
need to minimize |Im (t — f)| = s — q. This occurs when q = a, which imphes 
X is parallel to y. But 

S - a = Sr + Se- IVr + ye\ > i^r " l^rl) + " IVel) (23) 

by the triangle inequality, and both terms on the right are positive by (22). 
Therefore, s — q is maximized, for given and a^, when x,y^ and are 
all parallel. That is, the gain is maximized when the two antennas are in 
a line-of-sight configuration. Thus our Ansatz (20) about the dual complex 
structures Zg = Xe + ii/e and Zr = Xr — iyr of the emitting and receiving 
sources leads to an obviously valid conclusion. (It would lead to the opposite 
conclusion had we set Zr = Xr + iyr, which would also introduce unwanted 
singularities into the wavelet W{zr — z^) since yr — ^y+-) 

Note that W is invariant under complex spacetime translations: 

z'^ = Ze + C, z'^ = Zr + C =^ z! = z[. — z'^ = Z , ( = + if] G C^, 

which gives the new source vectors 

y'e = ye + ri, y'^ = yr- r], y' = + y'^ = y. (24) 

(The imaginary translation vector r] must be sufficiently small that y'^ and 
y'j. remain in the future cone.) Invariance under real spacetime translation 
is nothing new, but its extension to imaginary translations is new and in- 
teresting. Equation (24) states that the emitter-receiver pairs {ze,Zr) and 
{z'f,,z'j.) form equivalent channels. This is far from obvious, though it makes 
intuitive sense. If we do not use a line-of-sight configuration for the emitter 
and receiver, we must increase the aperture sizes Og = \yg\ and Ur = {y^l to 
obtain the same transmission quality. In particular, consider the following 
two "extreme" channels related to the channel depicted in Figure 1, which 
we will call channel A: 

1. Channel B, equivalent to A but with an idealized event receiver: 

C = iyr y'e = ye + yr, y'r = 0. 

2. Channel C, equivalent to A and B but with an idealized event emitter: 

C = -iye ^ ye = 0, y'r = ye + yr- 
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Figure 2: Channel B (top) has an ideahzed event receiver, while channel C 
(bottom) has an idealized event emitter. Both are equivalent to channel A 
of Figure 1. Since point emitters and receivers are omnidirectional, channel 
B needs a larger emitting antenna than channel A to make up for the lost 
directivity of the receiver, and channel C needs a larger receiving antenna 
than channel A to make up for the lost directivity of the emitter. 
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6 Bandwidths 



Even though our analysis is in the time domain, it is possible to assign 
handwidths to the emitting antenna, the receiving antenna, and the com- 
munication channel as a whole. We define these to be the reciprocals of the 
shortest pulse that can be emitted, received and exchanged, respectively. 
The bandwidth of the emitting and receiving antennas are defined as 



re(0) Se-fle' ' ^^(0) - Or ' 

and the effective channel bandwidth is, by (23), 

where a = \y\ = IVr + vd < + Cg is the effective channel aperture. 
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